Phase-tunable colossal magneto-heat resistance in ferromagnetic Josephson thermal valves 
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We propose a heat valve based on the interplay between thermal transport and proximity-induced exchange 
splitting in Josephson tunnel junctions. We demonstrate that the junction heat conductance strongly depends 
on the relative alignment of the exchange fields induced in the superconductors. Colossal magneto-heat resis- 
tance ratios as large as ~ 10 7 % are predicted to occur under proper temperature and phase conditions, as well 
as suitable ferromagnet- superconductor combinations. Moreover, the quantum phase tailoring, intrinsic to the 
Josephson coupling, offers an additional degree of freedom for the control of the heat conductance. Our pre- 
dictions for the phase-coherent and spin-dependent tuning of the thermal flux can provide a useful tool for heat 
management at the nanoscale. 
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The study of heat transport and dynamics in meso- [1] and 
nanoscopic [2] solid-state systems, is a research field that has 
attracted much attention in recent years because of the im- 
pressive progress achieved in nanoscience and nanofabrica- 
tion techniques. At such scale heat may play a significant role 
in determining the properties of the devices, and therefore it 
is of particular interest to control and manipulate [3, 4] the 
thermal flux as well to understand the origin of dissipative 
phenomena. Prototypical cases in which the understanding of 
heat transport is crucial are, for instance, the fine temperature 
control in ultrasensitive cryogenic radiation detectors [1], gen- 
eral cooling applications at the nanoscale [1], and the emerg- 
ing field of coherent caloritronic circuitry where the quantum 
phase allows for enhanced operation [5-11]. 

It has been known for a few decades that phase-dependent 
thermal transport through weakly-coupled superconducting 
condensates is in principle possible [12-16]. However, only 
recently the first Josephson heat interferometer was demon- 
strated [17-19]. The experiment of Ref. [19] proves that, 
in addition to the Josephson charge supercurrent, phase co- 
herence extends to dissipative observables such as the ther- 
mal current. This heat interferometer represents a prototypical 
building block to implement future coherent caloritronic cir- 
cuits like, for instance, heat transistors and thermal splitters. 

In this Letter we put forward the concept of a ferromag- 
netic Josephson junction acting as a thermal valve. In par- 
ticular, we address the interplay between thermal transport 
and proximity-induced exchange splitting in a Josephson tun- 
nel weak-link consisting of two superconducting electrodes 
with an internal exchange splitting. The latter is induced from 
nearby-contacted ferromagnetic layers [see Fig. 1(a)]. We 
show that the junction thermal conductance strongly depends 
on the relative alignment of the exchange fields induced in 
the superconductors. As a results, colossal magneto-heat re- 
sistance ratios as large as ~ 10 7 % are predicted to occur for 
suitable exchange fields and proper temperature conditions. 
Moreover, the quantum phase tailoring, characteristic for the 
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FIG. 1. (Color online) (a) A schematical view of the FSISF Joseph- 
son heat valve discussed in the text, (b) The exchange fields Qil,r) 
in the F layers are confined to the z — y plane, and are misaligned 
by an angle a. (c) Temperature dependence of the self-consistently 
calculated superconducting order parameter A for different values of 
the exchange field h. Aq is the zero-temperature, zero-exchange field 
order parameter and T c is the superconducting critical temperature. 



Josephson effect, adds a further degree of freedom for en- 
hanced heat conductance control. 

Our system is schematized in Fig. 1(a). It consists of two 
equal ferromagnet- superconductor bilayers (FSl,r) tunnel- 
coupled through an insulating barrier (I) and implementing a 
Josephson junction. The FSl and FSr bilayers are in thermal 
steady- state and reside at different temperatures Tl and Tr, re- 
spectively. For definiteness, we assume Tl > Tr so that the 
structure is temperature-biased only, while there is no volt- 
age drop across the Josephson junction, ts (tp) labels the S 
(F) layer thickness while (p denotes the macroscopic quantum 
phase difference over the junction. Furthermore, the z-axis is 
the one parallel to the magnetization (exchange field) of the 
left F layer Qil), which is kept fixed, whereas the one in the 
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right ferromagnet (Jir) is misaligned by an angle a [see Fig. 
1(b)]. Experimentally this can be achieved either by using fer- 
romagnetic films with different coercive fields or by pinning 
the magnetization in the left electrode through an exchange- 
bias with an additional magnetic layer [20]. Hr can therefore 
be freely rotated by applying an in-plane magnetic field as low 
as a few tens of Oe. 

We first derive an expression for the heat current (Q) flow- 
ing through the Josephson junction. The latter can be ex- 
pressed in terms of the quasiclassical Green's functions (GFs) 
of the left and right electrodes 
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Here the trace is taken over the spin(g)particle-hole space 
while the functions G R n) are 8x8 matrices in the 
Keldy sh(g)particle-hole(g) spin space : 
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where j = R,L. The symbols T and ? denote 4x4 matrices in 
particle-hole0spin and 2x2 matrices in spin-space, respec- 
tively. Furthermore, Rn is the normal-state resistance of the 
junction and e is the electron charge. 

We assume that the electrodes are in thermal equilibrium, 
thus the Keldysh component of the GFs is given by gf(s) = 
(gj —gj)Fj, where Fj = tanh[e/(27))] is the electronic distri- 
bution function, and 7) is the temperature of the j electrode. 
According to Eq. (1) there is a finite heat current flowing 
through the junction if Tr ^ Tl which is given by 



Q = f de.sTr [N L N R -M l M r cos <p] [F R -F L ] . 



(3) 



The two contributions to the heat current stem from the 
normal, Nj = (gj — gj )/2, and phase-coherent (anomalous), 
Mj = (fj — ff)/2, parts of the quasiparticle spectral function 
[12, 15]. Equation (3) is the generalization of the Maki- Griffin 
heat current equation [12] for the case of spin-dependent den- 
sity of states (DoS). In particular, we obtain the oscillatory 
behavior of the heat current as a function of the superconduct- 
ing phase difference (p predicted for the first time in Ref. [12], 
and recently demonstrated in Ref. [19]. We stress that a pure 
temperature bias across the junction is a crucial condition to 
preserve phase dependence in thermal transport. Indeed, any 
voltage drop occurring across the Josephson weak-link would 
make (p time-dependent and, therefore, the (p-dependent com- 
ponent of Q in Eq. (3) would not contribute to the DC heat 
transport [12, 14, 19]. 

Instead of analyzing the heat current, that depends on a 
generic temperature difference across the junction, we shall 
focus on the behavior of the thermal conductance (k) which is 
defined for small temperature differences as 
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FIG. 2. (Color online) (a) Magneto-heat resistance ratio MHR vs 
temperature T calculated for a few values of the exchange field h 
at (p = 0. (b) MHR ratio vs T calculated for the same values of 
h as in panel (a) at (p = n/2. (c) MHR ratio vs (p calculated for 
several values of the exchange field at T = 0AT c . (d) MHR ratio vs 
(p calculated at T = 0.5T C for the same values of h as in panel (c). 



where ST = T L -T R , and (dF/dT) = -e/[2T 2 cosh 2 (e/2T)}. 
By deriving the second equality we have assumed that ST <C 
T = (Tr + Tl)/2. Equations (1) and (4) are rather general, and 
allow to compute the heat current and the thermal conductance 
for an arbitrary tunneling junction provided that values of the 
GFs on both side of the interfaces are known. 

With the help of Eq. (4) we can determine the heat con- 
ductance for the junction sketched in Fig. 1(a). We assume 
that \hi\ = \h R \ =h, and that the S/F interface is highly trans- 
missive so that both the superconductor and the ferromag- 
net are strongly affected by the proximity effect [21, 22]. At 
the same time, in order to preserve superconductivity in the 
leads, we assume that the F layers are thin enough. In par- 
ticular, if the thickness tgrp) °f the superconducting (ferro- 
magnetic) layer is smaller than the characteristic length over 
which the GFs vary, one can integrate the quasiclassical equa- 
tions over the thickness of the S/F bilayers [24]. After such 
procedure one obtains for the retarded and advanced GFs 
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are the Pauli matrices in particle-hole space. We focus first on 
the case that the magnetizations of the F layers in Fig. 1 are 
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either parallel (P) or antiparallel (AP) to each other. Thus GFs 
g R ( A ) an( i f R ( A ) are 2 x 2 diagonal matrices in spin space with 
diagonal elements given by [24] 



fi = 
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where h and A are the effective values of the exchange 
field and superconducting order parameter in the S/F bilayer, 
respectively. In particular, A has to be determined self- 
consistently. The temperature dependence of the order pa- 
rameter for different values of h is shown in Fig. 1(c). The 
parameter T in Eqs. (5-6) accounts for the inelastic scattering 
energy rate within the relaxation time approximation [25-27] . 
Similar expressions hold for the advanced GFs by replacing 
in Eqs. (5-6) zT by — zT. The real part of the functions g± 
gives the modified DoS in the superconductors which is spin- 
dependent due to the finite exchange field in the F layers. 
The heat conductance is thus obtained from Eq. (4) 
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where X = P,AP, A P (e) = £ a=± [N L aN Ra - M La M Ra cos <p] 
and A AP (e) = 2 [N L +N R _ - M L+ M R _ cos <p] . We propose 
an experiment in which one can switch between the P and 
AP configurations, and determine the magneto-heat resistance 
(MHR) ratio defined as 



MHR 



KP ~ *AP 
K A P 
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In Fig. 2 we show the behavior of the MHR as a func- 
tion of temperature and the superconducting phase difference. 
All panels show an overall huge MHR ratio (~ 10 5 — 10 7 %) 
within a broad range of parameters. We demonstrate in this 
way that by switching between the P and AP configuration 
one realizes an almost perfect heat valve effect as the thermal 
conductance in the AP configuration is practically negligible 
with respect to that in the P one. This colossal MHR is one 
of the key results of the present letter. Figures 2(a) and 2(b) 
show that the heat valve effect is maximized at certain finite 
temperature (i.e., for T/T c ~ 0.1) and for sufficiently large 
exchange fields. Here T c is the superconducting critical tem- 
perature. It is worth emphasizing that due to the cos (p inter- 
ference term in Eq. (4) the MHR ratio can be additionally 
largely tuned by the phase difference between the supercon- 
ductors. Such a phase-tunable thermal transport mechanism 
originates from the Josephson effect and is unique to weakly- 
coupled superconductors [12]. In the lower panels of Fig. 2 
the MHR dependence on (p is displayed. The minimum value 
of the MHR is achieved for zero phase difference, whereas it 
reaches its maximum value for (p = 71. We also emphasize 
that the phase-coherent term in Eq. (4) does not describe pure 
tunneling of Cooper pairs [12, 13]. Furthermore, we point out 
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FIG. 3. (Color online) (a) MHR a ratio vs T calculated for several 
values of the misalignment angle a at h = 0.1 Aq and (p = 0. (b) 
MHR a ratio vs T calculated for the same a values as in panel (a) at 
h = 0.5Aq and (p = 0. (c) MHR ratio vs T calculated for a few values 
of r at h = 0.5A and (p = 0. (d) MHR ratio vs (p calculated fot the 
same T values as in panel (c) at h = 0.5Aq and T = 0.1T C . 



that while the P configuration maximizes the heat current, the 
DC Josephson effect is maximized by the AP one [24] . 

The obtained colossal MHR ratio can be understood by in- 
spection of Eq. (7). If we assume for simplicity that (p = 71 / 2 
[see Fig. 1(b)] then only the normal GFs contribute to Ax 
[cf. Eq. (4)]. The heat current (and hence the thermal con- 
ductance) is a consequence of quasiparticle transmission from 
the hot to the cold electrode. For a given energy, the num- 
ber of states available for the heat transport is given by the 
spectral function Ax which is the product of the DoS on both 
side of the tunnel barrier. Due to the exchange field, the DoS 
is spin-dependent, and shows a BCS-like shape with spin- 
dependent energy gap at A± = A ± h [see Eq. (5)], equivalent 
to a Zeeman-split superconductor in a magnetic field [28]. In 
the P configuration, the DoS of the left and right electrode 
coincide for both spin-up and spin-down, and therefore quasi- 
particles with energies around £ ~ A± contribute at most to 
the heat conductance [29]. The situation is different in the AP 
configuration, where the DoS for each spin-channel is shifted 
on both side of the barrier by an amount 2h. The main contri- 
bution to Kap comes from quasiparticles with energies £ ~ A + 
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and Aap is approximately a factor ~ yT/A smaller than yip. 
Moreover, in both the P and AP cases the contribution from 
Ap(ap) is weighted by the function edF /dT(e). The latter 
decays as e~ £ / 2T for £ > 2T and hence the main contribution 
to Kap in Eq. (7) (from £ ~ A + ) has an additional exponen- 
tially small factor e -/ */ r with respect to the main contribution 
to Kp (from £ ~ A_). All of this explains the smallness of k*ap 
and the huge MHR ratio obtained for sufficiently large values 
of the exchange field. 

As discussed above, the maximum MHR ratio is reached 
for a certain finite temperature. According to Figs. 2(a,b) a 
further increase of T leads to a decrease of the MHR, which 
can be explained, on the one hand, by the suppression of 
the energy gap A(T) and on the other hand, by the fact that 
increasing T the contribution from quasiparticles with ener- 
gies larger than A± becomes more and more important lead- 
ing to a smaller difference between Kap and Kp. Notice that 
for < (p < 71 / 2 the condensate part of the spectral func- 
tion Ax(£) [see Eq. (7)] gives a negative contribution to the 
heat conductance. This explains the lower values of MHR for 
small phase difference shown in Figs. 2(c) and 2(d). 

For an arbitrary angle a between the magnetizations of 
the left and right electrode [see Fig. 1(b)] we define MHR a 
as MHR a = (K a — K A p)/ K A p, where K a = k p cos 2 (oc/2) + 
Kap sin 2 ( a /2). The last expression for K a can be obtained 
straightforwardly from Eq. (1) by rotating the right Green 
function according to Gr = R a Go$ a , where Go is the GF for 
the case a = and_? a = exp[/T30i«/2] [23, 30]. In Figs. 3(a) 
and 3(b) we show the temperature dependence of MHR a for 
different values of a at <p = 0. All curves show similar behav- 
ior, and again very large values for the MHR can be achieved 
with a proper choice of the parameters. According to Figs. 
3(a) and 3(b), the effect is maximized for a = 0, i.e., when 
the junction is switched between the P and AP configurations. 
Figures 3(c) and 3(d) show the impact of the inelastic parame- 
ter r on the MHR. The overall tendency is that by increasing F 
the MHR ratio is reduced, as the "normal" character of trans- 
port is strengthened in the heat valve leading to a suppression 
of the large MHR ratio. The latter, indeed, originates from the 
presence of the superconducting gap. Moreover, as displayed 
in Fig. 3(c), the MHR ratio reaches its maximum at higher 
temperature by increasing T. 

In light of a realistic implementation of the present heat 
valve, soft ferromagnetic alloys such as Cui_ x Ni x [31] or 
Pdi_ x Ni x [32], which allow fine tuning of the exchange field 
through a suitable choice of x, combined with a conventional 
superconductor (e.g. aluminum or niobium) might be suit- 
able candidates. We note that all the results presented above 
have been obtained assuming highly-transparent S/F inter- 
faces. But nevertheless, they are qualitatively valid as well 
in the case of a finite S/F interface resistance R^. In such a 
case the superconductor still exhibits a spin- split DoS, how- 
ever with an additional damping factor determined by R^. The 
latter will suppress the MHR similarly as it does a finite T. 
Furthermore, according to our model one can also design the 



heat valve of Fig. 1(a) by using ferromagnetic insulators (FIs), 
as for example Eu chalcogenides barriers, instead of metallic 
ferromagnets. In such a case it was experimentally proved 
[33-35] that the DoS in the superconductor is modified and 
shows the spin- splitting needed to obtain the heat valve effect. 
Therefore, all the conclusions drawn above remain valid if one 
designs the junction by exploiting FIs for the F layers. 

With regards to potential applications, the present thermal 
valve can be used whenever a precise control and mastering 
of the temperature is required, for instance, for on-chip heat 
management as a switchable heat sink. This setup can be use- 
ful as well to tune the operation temperature of sensitive ra- 
diation detectors [1, 36]. In the context of quantum comput- 
ing architectures [37] the Josephson thermal valve can also be 
used to influence the behavior and the dynamics of two-level 
quantum systems through temperature manipulation. Simi- 
larly, the relation between Josephson critical supercurrent and 
the temperature can be exploited for designing tunable ther- 
mal Josephson weak-links of different kinds [1, 38, 39]. 

In conclusion, we have investigated thermal transport 
through a heat valve consisting of a Josephson junction be- 
tween two S/F bilayers as electrodes. In particular, we predict 
that the heat conductance depends strongly on the the rela- 
tive alignment of the magnetizations of the F layers. Under 
specific conditions of temperature bias and phase difference 
across the junction one can obtain a colossal magneto-heat 
resistance ratio as high as several orders of magnitude. The 
spin-dependent and phase-tunable mechanisms of heat flux 
control discussed in this letter will likely prove useful for ther- 
mal management at the nanoscale, and for the development of 
coherent spin caloritronic nanocircuits [40, 41] . 
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